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. We obtain in closed form the probability distribution for individual measurements of the stress 
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I energy tensor of two-dimensional conformal field theory in the vacuum state, smeared in time 

against a Gaussian test function. The result is a shifted Gamma distribution with the shift given 
by the previously known optimal quantum inequality bound. For small values of the central charge 
it is overwhelmingly likely that individual measurements of the sampled energy density in the 
vacuum give negative results. For the case of a single massless scalar field, the probability of 
finding a negative value is 84%. We also report on computations for four-dimensional massless 
scalar fields showing that the probability distribution of the smeared square field is also a shifted 
Gamma distribution, but that the distribution of the energy density is not. 
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The expectation value of the stress tensor of a quantum field, when suitably renormalized, 
can be used to give a description of the gravitational effects of quantum fields. However, this 
semiclassical theory clearly has limitations in that it does not describe fluctuations around 
the mean value. In recent years, several authors have discussed possible physical effects 
of these fluctuations. For recent reviews with further references, see Refs. |2|. Much of 
the work on this topic has involved seeking observables which can be related to integrals 
of the stress tensor correlation function. However, this amounts to calculating a variance, 
or a second moment, of the probability distribution for measurement results of the smeared 
stress-energy tensor. It is desirable to go further and determine the complete distribution 
function. In this paper, we will report results for such a function. 

Here we will be concerned with Minkowski spacetime in the vacuum state and seek the 
probability distribution function for a stress tensor operator which has been averaged with a 
sampling function. This averaging is crucial to having a well-defined probability distribution. 
(A recent attempt has been made to assign a distribution to a local, unsmeared, stress 
tensor operator |3|. However, in our view this is not a meaningful quantity.) There is 
a deep connection between the probability distribution for vacuum fluctuations, and the 
quantum inequalities (QIs) which place lower bounds on the expectation values of smeared 
stress tensors in arbitrary physically reasonable quantum states p-lO]. These lower bounds 



represent the lowest eigenvalues of the smeared operators. As a consequence of the Reeh- 
Schlieder theorem (see, e.g., they are also the lower bounds of the support of the 



corresponding probability distributions in the vacuum state [18[. In this paper, we will 
concentrate on conformal field theory (CFT) in two spacetime dimensions, but will also 
report a result for the square of a massless scalar field in four spacetime dimensions. Further 
details will be given in a later paper. 

We consider a general unitary, positive energy CFT with central charge c in two- 
dimensional Minkowski space (with H — signature and inertial coordinates x"). A general 



reference is 12j. Our notation is as follows: Tl and Tr are the left- and right-moving 



components of the (renormalized) stress-energy tensor, so that 



Too(a;°,xi) = -Tn(x°,xi) = Tfi(n)+Ti(i; 
Toi(x°,a;i)=Tio(x°,xi) = -Tn{u)+Ti 



\{v) (1) 



where u = — , v = + . As Tr and commute, we can find the probabil- 
ity distribution for a smearing of Tqo from the probability distributions for correspond- 
ing smear ings of Tl and Tr. We begin by considering Tr only and dropping the sub- 
script (identical results hold for Tl). The vacuum n-point functions will be denoted 
G„(Mn, . . . , Ml) = {T{un) ■ ■ -Tiuij) , with the convention that Go = (1) = 1; where one 
or more of the Uj is replaced by the letter /, this indicates smearing against a test function 
f{uj), and we write Qn[f] = Gn{f, •••,/) = {T{fY) for the n-th moment of T(/) in the 
vacuum state. All test functions will be assumed real- valued and rapidly decaying at infinity. 
Our goal is to determine an expression for the moment generating function 

MW] = E (2) 

n=0 

from which the underlying probability distribution can be determined under suitable con- 
ditions. Here, we take j f{u) du to be dimensionless and /i with dimensions of length^. We 
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start with the Ward identities in CFT (see, e.g., Sect. 3 of [l2|), which entail the recurrence 
relation 



G„(m„, ...,ui] 



n-l 

E 

j = l L 



C Gn^2{Un^l, ■ ■ ■ ,Uj, ■ ■ ■ ,Ui) 



{Un - Uj - iO)'^ 



2TT{Un — Uj — iO) TT{Un — Uj — zO)^ 



(3) 



where the hat denotes an omitted variable. Since Go = 1 and Gi{ui) = (T(mi)) = 0, it 
follows immediately that 

c 

G2iU2,Ul) = —^ 7^. (4) 

Smearing Eq. ([3]) against n copies of / and noting that the first line of the right-hand side 
provides n — 1 identical contributions, we have 



^„[/] = (n-l)^2[/]^n-2[/] + $^/n,,, 



where it may be shown that 



(5) 



(6) 



with the underbrace indicating that /* / is inserted in the j'th position from the right, and 
/ * / given after some computation by 



■oo 



27c(w — u) 



(7) 



Note that the integrand is nonsingular at w = u. In deriving this formula several steps of 
integration by parts have been used; no boundary terms are introduced provided / decays 
sufficiently fast at infinity. 

Now suppose that a one-parameter family fx{u) (A G M) can be found, satisfying the 
differential equation 

t^/.*/. (8) 

and initial condition /o = / (A has dimensions of length^). Then, replacing / by fx in 
Eq. ([5]), the sum over j can be rewritten as a derivative with respect to A [loj, yielding 



Gnifx] = in- l)g2[fx]Gn-2[fx] + ^Sn-l[fx]- 



(9) 



Multiplying by /i" ^/{n — 1)! and summing over n > 2, this becomes a partial differential 
equation 

gj^M[^ifx] = fig2[fx]M[^ifx] + g^M[f^fx] (10) 
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for M[fifx], on using the facts Qoifx] = 1, ^i[/a] = 0, and Gn[fJ'fx] = f^"'Qn[f\]- To solve this, 
we write M[f] = exp(W[f]), whereupon Eq. ( ITO|) becomes 



—w[fifx] = fig2[fx] + g^w[fih]. 



(11) 



As iy[/i/A] = for /i = 0, the unique solution is 



and we set A = to obtain 



dx (/i-A)e;2[/, 



A • 



(12) 



(13) 



Equation ([T3|) is similar to a formula derived in the Euclidean holomorphic picture (in 
contrast to our Lorentzian treatment) by Haba in Eq. (9) of [l3|. 

The differential equation Eq. ([8]) can be solved explicitly in the case of a Gaussian sam- 
pling function f{u) = e~'"^/'^^/(r-\/vr). Making the ansatz f\{u) = A{X)f{u), the calculation 



)(") = = ^M")- 



(14) 



reduces Eq. ([8]) to dA/dX = A(A)^/(r^7r), whose unique solution with ^4(0) = 1 is 



^(A) 



ITT 



-kt"^ — X 

We now have ^2 [/a] = A{Xfg2[f], so 

wM = g2[f] rdX{f,-x)A{xY 

Jo 

= G2[f] 



(15) 



r'^TT^ log 


( ^^^^ ) 






\7rr^ — jj J 





where we have used 



to compute 



24 
1 



log 



TTT 



(V — le] 



6 Jo 



G2[f] 



487r2 



247r2r4 ' 

The moment generating function for the Gaussian-averaged right-moving flux is then 

c/24 



(16) 

(17) 
(18) 



M[/i/] = e 



W[/if] 



(19) 
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It remains to determine the corresponding probability distribution P{uj) such that 



oo 



M[fif] = / dujP{u)e^''. (20) 



— oo 



As the moments obey the uniqueness conditions in the Hamburger moment theorem [14 
there can be only one such distribution. (One needs |a„| < CD^n\ for the n'th moment 
a„ = (i"M[/i/]/(i/i" 1^=0; which is satisfied for our case.) It is a shifted Gamma distribution, 
of form 

P{uj) = ^{u + a;o)^!i^i±^^^ exp(-/3(a; + c^o)), (21) 

r a 



with parameters 



and where ^ is the Heaviside function. Note that Eq. ( 1211) is singular at its lower limit 
provided that c < 24, but the singularity is integrable, so that the integrated probability is 
unity, as required. 

As described above, the lower bound of this distribution is also the best possible quantum 
inequality bound on the expectation value in a normalized state ip, so 

(T(/))^ > -coo (23) 

for all physically reasonable ip. This can be checked against the optimal bound given in 
Theorem 4.1 of |lo|. 

Thus we have the expected agreement between the two methods. The nonzero negative 
lower bound should be a general feature of stress tensor probability distributions, and is 
required by the vanishing mean value in the vacuum state, in contrast to the result claimed 
in Ref. Q. 

The cumulative distribution function of our probability distribution is an incomplete 
Gamma function. Perhaps surprisingly, it is overwhelmingly likely that an individual mea- 
surement will yield a negative result: the probability is 0.89 for c = 1 (independent of tq 
and r) and decreases as c increases, tending to 0.5 as c — )■ cxd. 

When both Tr and Ti are combined to obtain the overall probability distribution for the 
energy density Too(t,0), the effect is to replace c/24 by c/12 in the probability distribution 
and definition of coq. The averaged energy density operator is p = Ttt{t, 0) e~* dt . 

Then x = pr"^ is distributed according to the shifted Gamma distribution with probability 
density function 

_-c/12/ _|_ ^c/12-l 

P{x) = ^{x + xo) ^j^^ e-(^+-) , (25) 

where — xq/t^ is the quantum inequality bound on expectation values of Tu in arbitrary 
quantum states. The probability distribution is plotted in Fig. [1] for the case of a free 
massless scalar field in two-dimensional Minkowski spacetime, for which c = 1. In this case, 
the probability of a negative result is 0.84, and the corresponding QI bound coincides with 
that of S. 
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FIG. 1: The probability distribution P{x) for the smeared energy density of a massless scalar field 
in two-dimensional Minkowski spacetime is plotted. Here x = pr"^, where p is the energy density 
operator averaged in time with a Gaussian function of width r. The lower limit of P{x) occurs at 
X = —xq = — l/(127r), illustrated by the vertical line. 



n 





1 


2 


3 


4 


5 


6 


7 


8 




1 





2 


48 


1740 


83904 


5051640 


364724928 


30707616912 



TABLE I: Computed moments Mn = (47rr)^"^„[L] for the Lorentzian averaged squared field 
operator. 



This is a remarkable result which sheds light on the nature of vacuum energy fluctuations, 
at least in two-dimensional spacetime. The zero average value of the smeared energy density 
in the vacuum state arises, in a sequence of measurements, from many negative values being 
balanced by rarer, but larger positive values. 

Although the probability distribution for stress tensor operators in four-dimensional 
spacetime is not yet known, we have been able to construct distributions for the square 
of a massless scalar field by empirical methods. This amounts to calculating several mo- 
ments, and then showing that they fit to a shifted Gamma distribution to high accuracy. 
Using the Lorentzian sampling function L{t) = r/(7r(t^ + r^)), we computed the first twenty 
moments of J ■.^'^■.{t, x)L{t) dt in the vacuum state, using symbolic integration in Maple. 
After suitable normalization, the computed moments are integers; the first eight of which 
are presented in Table [H Remarkably, all twenty moments are fitted exactly by the shifted 
Gamma distribution Eq. ( 12T]) with parameters 



a 



1 

72' 



/3 



2^2 



967r2r2' 



(26) 



Because the form in Eq. (1211) is already normalized, we need the n = 1,2, and 3 moments 
to fit the above parameters. However, as the remaining 17 moments are reproduced exactly, 
this provides strong evidence that this is the correct distribution. 

The lower bound of the probability density support is —ujq, which provides a conjectured 
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optimal QI bound 

I I O)^)L(t) dt > (27) 

for all normalized ijj in the relevant domain. By comparison, the methods of yield a 

QI H 



{iP I ■.^'■.{t,x)ilj)f{t)dt > 



r2 



d 



dt^^ 



dt. (28) 



Substituting the Lorentzian L in place of /, the lower bound is — l/(64'7r^r^), which is weaker 
by a factor of 3/2 than (and therefore does not contradict) the conjectured bound just given. 

Just as in the case of fluxes and energy densities in two dimensions, the probability 
distribution diverges integrably at its lower limit. In this case, the probability of a negative 
outcome for a measurement of the averaged squared field is 0.95. 

The first fifteen moments of :$^: smeared against the squared Lorentzian S{t) = 
2r'^/(7r(t^ + T^)^) have also been computed and are exactly fitted by the shifted Gamma 
distribution f l?T]) with parameters a = 1/45, (3 = 87r^r^/15, and Uq = l/(247r^r^). The cor- 
responding conjectured optimal QI bound is again stronger by a factor of 3/2 than that of 
Eq. ( 128|) : it is tempting to conjecture that the optimal bound for general sampling functions 
is given by ( |28|) with a 12 replacing the 8. 

As in the case of the stress tensor in CFT with c < 24, these distributions for the squared 
scalar field are singular at the lower limit, although with an integrable singularity. The 
physical effect of this singularity is to favor negative fiuctuations close to the lower bound. 
Note that there nothing unphysical about this singularity, as the observable quantity, the 
probability of an outcome in a finite interval, is finite. 

We have also attempted a similar fit to the Lorentzian average of the squared time 
derivative of a massless scalar field and to the squared electric field. In both shifted 
Gamma distribution which is fit to the lower moments underestimates the higher moments. 
This suggests that the correct probability distribution for these quantities, and for the four- 
dimensional stress tensor, has a positive tail which falls off more slowly than is described by 
the shifted Gamma distribution. 

The probability distribution for the energy density in four- dimensional theories is of con- 
siderable interest. One application is to inflationary cosmology, where quantum stress tensor 
fluctuations might contribute a potentially observable component to the cosmological den- 



sity fluctuations [15[. This component would be non-Gaussian in a way which is associated 
with the skewness of the quantum stress tensor probability distribution. Another, perhaps 
more exotic, application to cosmology arises in models which employ anthropic reasoning to 
attempt to compute probabilities of various observables. These models require a counting 
of observers, possibly including "Boltzmann brains" which have nucleated from the vacuum 
in deSitter or even Minkowski spacetime. If this is the more prevalent type of observer, 
it would greatly complicate atternpts at anthropic prediction. (For further discussion and 
references, see, for example Refs. |16l. Il7l|.) The key to studying this question is in the details 
of the long positive tail of the probability distribution. 
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[18] The connection begins with the observation that the following three statements are equivalent, 
where A is a self-adjoint operator (with domain D{A) and spectrum cr{A)) representing a 
quantum observable: (a) {tp \ Ail^) > —Q for all normalised tp G F){A); (b) cr{A) C [—Q,oo); 
(c) the probability measure //^ that describes the distribution of individual measurement 
results of A in state ip has support supp/x^ C [—Q,oo). [The implication (a) =^ (b) is the 
Rayleigh-Ritz variational principle; (b) =^ (c) because measurements always give values in 
the spectrum; (c) =^ (a) because any random variable taking values in [—Q, oo) has its 
expected value in that range.] In general, supp fi^p (for a fixed choice of ^p) can be a proper 
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subset of (^{A), and in particular the infimum of the support of can exceed infer (74). For 
example, this happens if inf cr(A) is an isolated eigenvalue whose eigenfunction is orthogonal 
to ip. However, the spectrum and the support of the probability measure coincide if ^ is a 
local observable, such as a smeared field, and "0 is the vacuum state of a Minkowski space 
QFT obeying the standard assumptions of algebraic quantum field theory 111]. To see this, 
note first that if A G (^{A) then P^, the spectral projection of A corresponding to the interval 
(A — e, A + e), is nonzero for any e > 0. As is a local observable (because A is) the Reeh- 
Schlieder theorem entails that it cannot annihilate the vacuum, and hence there is a nonzero 
probability of obtaining a measurement of A in the interval (A — e, A + e), which must 

therefore intersect the support of fiQ. Taking e — )• we conclude that A G supp/U^ (which is a 
closed subset of M). Varying A, this argument establishes that cr{A) C supp/i^, and we already 
know that the reverse inclusion holds [by the implication (b) =^ (c) above]. In particular, 
cr(A) and supp/Xf^ have a common infimum, which provides a sharp quantum inequality bound 
(■0 ] Alp) > inf supp/if7 for all normalized ip £ D{A) by our earlier observations. 

[19] For example, use /a+<5 = fx + Sfx * fx + OiS"^) and linearity of in each slot. 

[20] Compare with the QI given in Eq. (5.5) of [9|; the starting point for Eq. ([28j) is to set p = \/y/oj 
in Eq. (3.11) of \ 
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